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ABSTRACT 



Aims. We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like 
oscillating stars to compare them with recent observations 

Methods. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to char- 
acterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuations in the 
uppermost part of the convective zone of such stars. These constraints, coupled with a theoretical model for stochastic 
excitation, provide the rate 7 at which energy is injected into the p-modes by turbulent convection. These energy rates 
are compared with those derived directly from the 3D simulations. 

Results. The excitation rates obtained from the 3D simulations are systematically lower than those computed from the 
semi-analytical excitation model. We find that f m „, the 7 maximum, scales as (L/M) s where s is the slope of the 
power law and L and M are the mass and luminosity of the ID stellar model built consistently with the associated 
3D simulation. The slope is found to depend significantly on the adopted form of Xh> the eddy time-correlation; using 
a Lorentzian, x\, results in s = 2.6, whereas a Gaussian, Xk j gives s = 3.1. 

Finally, values of U max , the maximum in the mode velocity, are estimated from the computed power laws for Tmax 
and we find that Knax increases as (L/M) sv . Comparisons with the currently available ground-based observations show 
that the computations assuming a Lorentzian \k yield a slope, sv, closer to the observed one than the slope obtained 
when assuming a Gaussian. We show that the spatial resolution of the 3D simulations must be high enough to obtain 
accurate computed energy rates. 

Key words, convection - turbulence - stars:oscillations - Sumoscillations 



"-j 1. Introduction 



Stars with masses M < 2M Q have upper convective zones 
where stochastic excitation of p-modes by turbulent con- 
vection takes place as in the case of the Sun. As such, these 
stars are often referred to as solar-like oscillating stars. One 
of the m ajor goals of the future space sei smology mission 
CoRoT (|Baglin fc The Corot Team! Il998l) . is to measure 
the amplitudes and the line-widths of these stochastically 
driven modes. From the measurements of the mode line- 
widths and amplitudes, it is possible to inf er the rates at 
which acoustic modes are excited (see e.g. [Baudin et all 
2005). Such measurements will then provide valuable con- 
straints on the theory of stellar oscillation excitation and 
damping. In turn, improved models of excitation and damp- 
ing will provide valuable information about convection in 
the outer layers of solar-like stars. 

The mechanism of stochastic excitation has been 
modeled by seve r al authors (e.g. iGoldreich fc Keelev , 
19771 lOsakil. 11990) : iBalmforthl . 119921 : IGoldreich et all 11 994: 
Samadi fc Goupil 12001. for a review see lStein et al. [ [2001 . 



are excited by turbulent convection but require an accu- 
rate knowledge of the time averaged and - above all - the 
dynamic properties of turbulent convection. 



Eddy 



time-correlations. In 



the 



These models yield the energy rate, CP, at which p-modes 
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approach of 

ISamadi fc Goupill (|200lL hereafter Paper I), the dy- 
namic properties of turbulent convection are represented 
by Xk, the frequency component of the auto-correlation 
product of the turbulent velocity field; Xfc c an be related 
to the convective eddy time-correlations. I Samadi et alj 
(|2003bl hereafter Paper III) have shown that the Gaussian 
function usually used for modeling Xk is inappropriate and 
is at the origin of the under-estimation of the computed 
maximum value of the solar p-modes excitation rates when 
compared with the observations. On the other hand, the 
authors have shown that a Lorentzian profile provides the 
best fit to the frequency dependency of Xk as inferred from 
a 3D simulation of the Sun. Indeed, values of CP computed 
with the model of stochastic excitation of Paper I and using 
a Lorentzian for Xk — Xk 1S better at reproducing the solar 
seismic observations whereas a Gaussian function, Xfe: 
under-estimates the amplitudes of solar p-modes. Provided 
that such a non-Gaussian model for Xk is assumed, the 
model of stochastic excitation is - for the Sun - rather 
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satisfactory. An open question, which we address in the 
present paper, is whether such non-Gaussian behavior 
also stands for other solar-like oscillating stars and what 
consequences arise for the theoretical excitation rates, CP. 



Stochastic excitation in stars more luminous than the Sun. 

In the last five years, solar- like oscillations have been 
det ected in several sta r s (se e for instance the review 
by iBedding fe Kieldsenl ([2003)). Theoretical calculations 
result in an overe s timat i on of their amp li tudes (see 
iKieldsen fc Bedding) . 120011: iHoudek fc Goughl . 12002ft . For 



instance, using Goug h's (|1976t I 1973) non-local and time 



dependent treatment of convection, IHoudek et alj ([1999) 
have calculated expected values of VJ nax , the maximum 
oscillation amplitudes, for different solar-like oscillating 
stars. Their calculations, based on a simplified excitation 
model, imply that V max of solar-type oscillations scale as 
(L/M) 1 - 5 where L and M are the lum inosity and mass 
of the star (see IHoudek fc Goughl . 120021 hereafter HG02 
). A similar scaling l aw w as empirically found earlier by 
IKieldsen fc Bedding! (|1995f ). As pointed out by HG02, all 
these scaling laws overestimate the observed amplitudes of 
solar-like oscillating stars hotter and more massive than 
the Sun (e.g. /3Hydri, r^Bootis, Procyon, £ Hydrae). As the 
mode amplitude results from a balance between excitation 
and damping, this overestimation of the mode amplitudes 
can be attributed either to an overestimation of the exci- 
tation rates or an underestimation of the damping rates. 
In turn, any overestimation of the excitation rates can be 
attributed either to the excitation model itself or to the 
underlying convection model. 

All the related physical processes are complex and diffi- 
cult to model. The present excitation model therefore uses 
a number of approximations such as the assumption of in- 
compressibility, and the scale length separation between the 
modes and the turbulent eddies exciting the modes. It has 
been shown that the current excitation model is valid in 
the case of the Sun (Paper III) , but its validity in a broader 
region of the HR-diagram has not been confirmed until now. 

Testing the validity of the theoretical model of stochas- 
tic excitation with the help of 3D simulations of the outer 
layers of stellar models is the main goal of the present paper. 
For that purpose, we compare the p-mode excitation rates 
for stars with different temperatures and luminosities as 
obtained by direct calculations and by the semi-analytical 
method as outlined below. 

Numerical 3D simulations enable one to compute di- 
rectly the excitation rates of p-modes for stars with 
various temperatures and luminosities. For instance this 
was a lready undertaken for the Sun by IStein fc Nordlundl 
(12001ft using the num erical approach introduced in 
iNordlund fc Steir] (|200lh . Such calculations will next be 
called "direct calculations" . They are time-consuming and 
do not easily allow massive computations of the excitation 
rates for stars with different temperatures and luminosities. 
On the other hand, an excitation model offers the advantage 
of testing separately several properties entering the excita- 
tion mechanism which are not well understood or modeled. 
Furthermore, once it is validated, it can be used for a large 
set of ID models of stars. 

As it was done for the Sun in ISamadi et all (|2003d 
hereafter Paper II) and Paper III, 3D simulations can also 
provide quantities which can be implemented in a formula- 



tion for the excitation rate CP, thus avoiding the use of the 
mixing-length approach with the related free parameters, 
and assumptions about the turbulent spectra. Such calcu- 
lations will next be called "semi-analytical calculations". 

We stress however that in any case, we cannot avoid the 
use of ID models for computing accurate eigen-frequencies 
for the whole observed frequency range. In the present pa- 
per, the ID models are constructed to be as consistent as 
possible with their corresponding 3D simulations, as de- 
scribed in Sect. 3. 

This paper is organized as follows: In Sect. 2 we 
present the methods considered here for computing CP, that 
is the so-cal led "direct" method based on Nordlund & 
Stein's (2001) approach (Sect. 2.1) and the so-called "semi- 
analytical" method based on the approach from Paper I, 
with modifications as presented in Papers II & III and in 
the present paper (Sect. 2.2). 

Comparisons between direct and semi-analytical calcu- 
lations of the excitation rates are performed in seven rep- 
resentative cases of solar-like oscillating stars. The seven 
3D simulations all have the same number of mesh points. 
Sect. 3 describes these simulations and their associated ID 
stellar models. 

The 3D simulations provide constraints on quantities re- 
lated to the convective fluctuations, in particular the eddy 
time-correlation function, Xk, which, as stressed above, 
plays an important role in the excitation of solar p-modes. 
The function Xk is therefore inferred from each simulation 
and compared with simple analytical function (Sect. 4). 

Computations of the excitation rates of their associ- 
ated p-modes are next undertaken in Sect. 5 using both 
the direct approach and the semi-analytical approach. In 
the semi-analytical method, we employ model parameters 
as derived from the 3D simulations in Sect. 4. 

In Sect. 5.2 we derive the expected scaling laws for CP max , 
the maximum in CP, as a function of L/M with both the di- 
rect and semi-analytical methods and compare the results. 
This allows us to investigate the implications of such power 
laws for the expected values of V max and to compare our re- 
sults with the seismic observations of solar-like oscillations 
in Sect. 5.3. We also compare wi t h previous theoretical re- 
sults (e.g. IKieldsen fc Bedding) . 119951 : IHoudek fc Goughl 
12002ft . 

We finally assess the validity of the present stochastic 
excitation model and discuss the importance of the choice 
of the model for \k m Sect. 6. 



2. Calculation of the p-mode excitation rates 

2.1. The direct method 

The energy input per unit time into a given stellar acoustic 
mode is calculate d num erically according to Eq. (74) of 
INordlund fc Steinl (|200lD multiplied by S, the area of the 
simulation box, to get the excitation rate (in erg s _1 ) : 



CP(w ) 



8 Av E u 



dr AP„ ad (r, w ) 



96 
dr 



(1) 



where AP nac i(r, ui) is the discrete Fourier component of the 
non-adiabatic pressure fluctuations, AP na( j(r, t), estimated 
at the mode eigenfrequency ujq = 2ttvq 1 £ r is the radial com- 
ponent of the mode displacement eigenfunction, Av = l/T 
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the frequency resolution corresponding to the total simula- 
tion time T a nd £^ n is the mode energ y per unit surface 
area denned in lNordlund k Steinl (|2001l their Eq. (63)) as: 



f — _ 



i"*/**'(b) ! 



(2) 



Note that Eq. ([T|) corresponds to the direct calculation of 
PdV work of the non-adiabatic gas and turbulent pressure 
(entropy and Reynolds stress) fluctuations on the modes. 
The energy in the denominator of Eq. (TT|) is essentially the 
mode mass. The additional factor which turns it into en- 
ergy is the mode squared amplitude which is arbitrary and 
cancels the mode squared amplitude in the numerator. For 
a given driving (i.e. PdV work), the variation of the mode 
energy is in yersely proportional to t he mode energy (see 
Sect. 3.2 of iNordlund k Steinl . |2001[ ). Hence, for a given 
driving, the larger the mode energy (i.e., the mode mass or 
mode inertia) the smaller the excitation rate. 

In Eq. (1) the non-adiabatic Lagrangian pressure 
fluctuation, AP na d(r, u>), is calculated as the following: 
We first compute the non-adiabatic pressure fluctuations 
AP na d(r, t) according to Eq. IA.3l in Appendix A. We then 
perform the temporal Fourier transform of AP na d (r, t) at 

each depth r to get AP na( j (r, uj) . 

The mode displacement eigenfunction £ r (r) and the 
mode eigenfrequency luq are calculated as explained in 
Section 3. Its vertical derivative, d^ T /dr, is normalized by 
the mode energy per unit surface area, £ Wo , and then multi- 
plied by AP na d. The result is integrated over the simulation 
depth, squared and divided by 8 Av. We next multiply the 
result by the area of the simulation box (§) to obtain f, the 
total excitation rates in erg s _1 for the entire star. Indeed 
the nonadiabatic pressure fluctuations are uncorrelated on 
large scales, so that average AP^ ad is inversely proportional 
to the area. Multiplication by the area of the stellar simu- 
lation gives the excitation rates for the entire star as long 
as the domain size is sufficiently large to include several 
granules. 

2.2. The semi-analytical method 

Calculations of excitation rates by the semi-analytical 
method are based on a model of stochastic excitation. The 
excitation model we consider is the same as presented in 
Paper I. In this mo del of excitation and i n con t rast to pre- 



vious models (e.g. Goldreich k Keelevl 119771 : iBalmforthl 



Il992t I Goldreich et all . 1 1994) . the driving by turbulent con 
vection is ensured not only by the Reynolds stress tensor 
but also by the advection of the turbulent fluctuations of 
entropy by the turbulent movements (the so-called entropy 
source term). 

As in Paper I, we consider only radial p-modes. We 
do not expect any significant differences for low £ degree 
modes. Indeed, in the region where the excitation takes 
place, the low £ degree modes have the same behavior as the 
radial modes. Only for very high t degree modes (£ 3> 100) 
- which will not be seen in stars other than the Sun - can a 
significant effect be expected, as is quantitatively confirmed 
(work in progress). 

The excitation rates are computed as in Paper II, except 
for the change detailed below. The rate at which a given 
mode with frequency ujq = lixv^ is excited is then calculated 
with the set of Eqs. (l)-(ll) of Paper II. These equations 



are based on the excitation model of Paper I, but assume 
that injection of acoustic energy into the modes is isotropic. 
However, Eq. (10) of Paper II must be corrected for an 
analytical error (see ISamadi et a"il . l2005t) . This yields the 
following correct expression for Eq. (10) of Paper II: 



Sr(t, u ) 



00 dk E(k,r) E{k,r) 



U 







II 







-t-oo 



(3) 



where uq — t/WJSu, $ is Gough's (|l977l ) anisotropy factor, 
u is the rms value of it, the turbulent velocity field, k the 
wavenumber and Xfc( w ) is t Qe frequency component of the 
correlation product of u. 

The method then requires the knowledge of a number 
of input parameters which are of three different types: 

1) Quantities which are related to the oscillation modes: 
the eigenfunctions (£ r ) and associated eigen- frequencies 

2) Quantities which are related to the spatial and time 
averaged properties of the medium: the mean density 
(po)i Q s = ((dp/ds) ) - where s is the entropy, p the 
gas pressure and (. . .) denotes horizontal and time av- 
erages - the mean square of the vertical component of 
the convective velocity, (w 2 ), the mean square of the 
entropy fluctuations, (s 2 ), and the mean anisotropy, $ 
(Eq. (2) of Paper II). 

3) Quantities which contain information about spatial and 
temporal auto-correlations of the convective fluctua- 
tions: the spatial spectrum of the turbulent kinetic en- 
ergy and entropy fluctuations, E(k) and P s (fc), respec- 
tively, as well as the temporal spectrum of the correla- 
tion product of the turbulent velocity field, Xfc- 

Eigen-frequencies and eigenfunctions [in 1) above] are 
computed with the adiabatic pulsat i on co de ADIPLS 
(jChristensen-Dalsgaard k Berthomieul . Il99ll ) for each of 
the ID models associated with the 3D simulations (see Sect. 

The spatial and time averaged quantities (in 2) and 3) 
above) are obtained from the 3D simulations in the manner 
of Paper II. For E(k), however, we use the actual spectrum 
as calculated from the 3D simulations and not an analyt- 
ical fit as was done in Paper II. However as in Paper II, 
we assume for E s (k) the fc-dependency of E(k) (we have 
checked this assumption for one simulation and found no 
significant change in J 1 ). 

For each simulation, we determine Xk as in Paper III (c/. 
Sect. 2]). Each Xk is then compared with various analytical 
forms, among which some were investigated in Paper III. 
Finally we select the analytical forms which are the closest 
to the behavior of Xk and use them, in Section^ to compute 

y. 

3. The convection simulations and their associated 
ID models 

Numerical simulations of surface conve ction for seven dif- 
ferent solar-like stars were performed bv lTrampedach et al.l 
(1999). These hydrodynamical simulations are character- 
ized by the effective temperature, T e ff and acceleration of 
gravity, g, as listed in Table[TJ The solar simulation with the 
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Fig. 1. Location of the convection simulations in the HR 
diagram. The symbol sizes vary proportionally to the stellar 
radii. Evolutionary tracks of stars, with masses as indicated, 
were calculated on the ba se of Christensen-Dals g aard's 
stellar evolutionary code (IChristensen-Dalsgaardl [1982; 
Christcnsen-Dalsgaard fc Frand scn. 1983^ 



also by using T-t r elation s derived from the simulations 
(|Trampedach et al.l . l2006rj) . 

The average stratifications of the 3D simulations, 
augmented by the fitted ID envelope models in the 
interior, were used as the basis for the cigcnmodc 
calculations using the adiabatic pulsati on code by 
IChristensen-D alsgaar cTfe Berthomieul (|l991). These com- 
binations of averaged 3D simulations and matched ID en- 
velope models will, from hereon, be referred to as the ID 
models. 

The positions of the models in the HR diagram are pre- 
sented in Figure [1] and their global parameters are listed 
in Table [U Five of the seven models correspond to ac- 
tual stars, while Star A and StarB are merely sets of at- 
mospheric parameters; their masses and luminosities are 
therefore assigned somewhat arbitrarily (the L/M-ratios, 
only depending on T c g and g, are of course not arbitrary). 



star 


Toff 
[K] 


M/M© 


R/Rq 


L/L Q 


LMq/MLq 


q Cen B 


5363 


0.90 


0.827 


0.51 


0.56 


Sun 


5802 


1.00 


1.000 


1.02 


1.02 


Star A 


4852 


0.60 


1.150 


0.66 


1.10 


q Cen A 


5768 


1.08 


1.228 


1.50 


1.38 


StarB 


6167 


1.24 


1.769 


4.07 


3.28 


Procyon 


6470 


1.75 


2.102 


6.96 


3.98 


r\ Boo 


6023 


1.63 


2.805 


9.31 


5.71 



Table 2. Fundamental parameters of the ID-models asso- 
ciated with the 3D simulations of Table Q] 



same input physics and number of mesh points is included 
for comparison purposes. The surface gravity is an input 
parameter, while the effective temperature is adjusted by 
changing the entropy of the inflowing gas at the bottom 
boundary. The simulations have 50x50x82 grid points. All 
of the models have solar-like chemical composition, with 
hydrogen abundance X = 0.703 and metal abundance 
Z = 0.0145. The simulation time-series all cover at least 
five periods of the primary p-modes (highest amplitude, 
one node at the bottom boundary), and as such should be 
sufficiently long. 

The convection simulations are shallow (only a few 
percent of the stellar radius) and therefore contain only 
few modes. To obtain mode eigenf unctions, the simu- 
lated domains are augmented by ID envelope models in 
the interior by means of the ste l lar en velope code by 
lChristensen-D alsgaar d~fc Frandsenl (|1983af ). Convection in 
the envelope models is ba sed on the mixing-length formal- 
is m (lBohm-Vitensd.[l958l). 

iTrampedach et ail (|2006af ) fit ID stellar envelopes to 
average stratifications of the seven convection simulations 
by matching temperature and density at a common pres- 
sure point near the bottom of the simulations. The fit- 
ting parameters are the mixing-length parameter, a, and 
a form-factor, /3, in the expression for turbulent pressure: 
^turb = Pq u muT' where umlt is the convective veloci- 
ties predicted by the mixing-length formulation. A con- 
sistent matching of the simulations and ID envelopes is 
achieved by using t he same equation of state (EOS) by 
iDappen et ail (| 19881 also referred to as the MHD EOS, with 
reference to Mihalas, Hummer, a nd Dappen) and opac- 
ity distribution functions (ODF) bv lKuruc 3 (ll992al lbT). and 



4. Inferred properties of Xk 

For each simulation, Xk{w) is computed over the whole 
wavenumber (fc) range covered by the simulations and at 
different layers within the region where modes are excited. 
We present the results at the layer where the excitation is 
maximum, i.e., where uq is maximum, and for two repre- 
sentative wavenumbers: fc = fc max at which E(k) peaks and 



10fc n 



where fcmin is the first non-zero wavenumber 



of the simulations. Indeed, the amount of acoustic energy 
going into a given mode is largest at this layer and at the 
wavenumber fc ~ fcmax, provided that the mode frequency 
satisfies: loq < (fcmax Uq)- Above loq ~ fcmax M , the effi- 
ciency of the excitation decreases rapidly. Therefore low 
and intermediate frequency modes (i.e., u>q < k max Uo) are 
predominantly excited at fc ~ k max . On the other hand, 
high frequency modes are predominantly excited by small- 
scale fluctuations, i.e. at large fc. The exact choice of the 
representative large wavenumber is quite arbitrary; how- 
ever it cannot be too large because of the limited number 
of mesh points fc < 25 fc m m and in any case, the excita- 
tion is negligible above fc ~ 20 fc m in- We thus chose the 
intermediate wavenumber fc = 10 fc m in- Fig. [2] presents Xk 
as obtained from the 3D simulations of Procyon, aCenB 
and the Sun, at the layer where uq is maximum and for 
the wavenumber fc max . Although defined as a function of u>, 
for convenience, \k is plotted as a function of v = lu/2tt 
throughout this paper. Fig. [3] displays Xk for fc = 10 fc m in- 
Results for the other simulations are not shown, as the re- 
sults for Procyon, a CenB and the Sun correspond to three 
representative cases. 
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star 


^sim 


size 


logs 




Hp 


Lh/Hp 


C s 


ts 


^sim / 




[min] 


[Mm 3 ] 




[K] 


[km] 




[km s" 1 ] 


N 




q Cen B 


59 


4.0x 4.0x 2.2 


4.5568 


5363 


95. 


42.1 


7.49 


12.72 


278.3 


Sun 


96 


6.0x 6. Ox 3.4 


4.4377 


5802 


134 


44.8 


7.78 


17.30 


332.9 


Star A 


80 


11. 6x11. 6x 6.4 


4.0946 


4851 


316 


36.7 


7.98 


39.66 


121.0 


a Cen A 


44 


8.9x 8.8x 5.1 


4.2946 


5768 


189 


47.1 


7.81 


24.17 


109.2 


StarB 


110 


20.7x20.7x11.3 


4.0350 


6167 


359 


57.7 


7.76 


46.29 


142.6 


Procyon 


119 


20.7x20.7x10.9 


4.0350 


6470 


380 


54.5 


7.52 


50.50 


141.4 


77 Boo 


141 


36.9x36.9x16.3 


3.7534 


6023 


709 


52.0 


7.40 


96.13 


88.0 



Table 1. Characteristics of the convection 3D simulations: t s - lm is the duration of the relaxed simulations used in the 
present analysis, H p is the pressure scale height at the surface, the size of the box in the horizontal direction, C s the 
sound speed and t s the sound travel time across H p . All the simulations have a spatial grid of 50x50x82. 



In practice, it is not easy to implement directly in the 
excitation model the v- variation of Xk inferred from the 3D 
simulations. An alternative and convenient way to compute 
3 is to use simple analytical functions for Xk which are cho- 
sen so as to best represent the 3D results. We then compare 
Xk computed with the 3D simulations with the following 
simple analytical forms: the Gaussian form 



1 



the Lorentzian form 
l, \ 1 1 

TTWfe/2 l + (2u/LU k f 

and the exponential form 



1 £ -\2u/u k \ 



(4) 



(5) 



(G) 



In Eqs.(|3][()]), u>k is the line- width of the analytical function 
and is related to the velocity u k of the eddy with wave 
number k as: 



LL>k — 2 kuk 



(7) 



In Eq. Q, Uk is calcu lated from the kinetic energy spectrum 
E(k) as lsteinl . [l96l 



_>/.- 



dkE(k) 



(8) 



As shown in Fig. [2] and Fig. O the Lorentzian x\ does 
not reproduce the ^-variation of Xk satisfactorily. This is 
particularly true for the solar case. This contrast with the 
results of Paper III where it was found that reproduces 
nicely - at the wavenumber where E is maximum - the 
^-variation of Xk inferred from the solar simulation inves- 
tigated in paper III. These differences in the results for 
the solar case can be explained by the low spatial resolu- 
tion of the present solar simulation compared with that of 
Paper III. Indeed we have compared different solar simu- 
lations with different spatial resolution and found that the 
^-variation of Xk converges to that of x\ as the spatial res- 
olution increases (not shown here). This dependency of Xk 
with spatial resolution of the simulation is likely to hold for 
the non-solar simulations as well. This result then suggests 
that Xk is in fact best represented by the Lorentzian form, 

As a consequence, realistic excitation rates evaluated di- 
rectly for a convection simulation should be based on sim- 
ulations with higher spatial resolution. However the main 



goal of the present work is to test the excitation model, 
which can be done with the present set of simulations. 
Indeed, we only need to use as inputs for the excitation 
model the quantities related to the turbulent convection 
(E(k), Xk,---) as they are in the simulations, no matter 
how the real properties of Xk are. 

For the present set of simulations, we compare three an- 
alytical forms of Xk- Lorentzian, Gaussian and exponential. 
For large k, Xk is overall best modeled by a Gaussian (see 
Fig.[3]for k = 10 fc m in). For small k (see Fig.[5]for k = fc max ) 
both the exponential and the Gaussian are closer to Xk than 
the Lorentzian. 

For a given simulation, depending on the frequency, dif- 
ferences between Xk(v) and the analytical forms are more 
or less pronounced. 

The discrepancy between Xk(y) inferred from the 3D 
simulations and the exponential or the Gaussian forms vary 
systematically with stellar parameters; decreasing as the 
convection gets more forceful, as measured by, e.g., the 
turbulent- to total-pressure ratio. Of the three simulations 
illustrated in Fig. [21 Procyon has the largest and a Cen B 
has the smallest Pturb/-ftot-ratio. 

As a whole for the different simulations and scale lengths 
k, we conclude that the ^-variation of Xk in the present set 
of simulations lies between that of a Gaussian and an ex- 
ponential. However, neither of them is c ompletely satisfac- 
tory. A ctually a recent detailed study bv lGeorgobiani et alj 
(2006], in preparation) tends to show that %fe cannot sys- 
tematically be represented at all wavenumbers by a simple 
form such as a Gaussian, an exponential or a Lorentzian, 
but rather needs a more generalized power law. Hence, more 
sophisticated fits closer to the simulated ^-variation of Xk 
could have been considered, but for the sake of simplicity 
we chose to limit ourselves to the three forms presented 
here. 



5. p-mode excitation rates across the HR diagram 

5.1. Excitation rate spectra (7(f)) 

For each simulation, the rates CP at which the p-modes of 
the associated ID models are excited are computed both di- 
rectly from the 3D simulations and with the semi-analytical 
method (see Sect. [2]). In this section, the semi-analytical 
calculations are based on two analytical forms of Xk'- a 
Gaussian and an exponential form as described in Sect. 4. 
The Lorentzian form as introduced in Sect.[4]is not investi- 
gated in the present section. Indeed our purpose here is to 
test the model of stochastic excitation by using constraints 



Samadi R. et al.: Excitation of solar-like oscillations across the HR diagram 



1 .000 




v = cj/2-n [mHz] 









Sun : 




Cj^,. . 




r 


V , 












— Gaussian 

Lorentzian 

Exponential 


#\ 

T\\ - 




v = co/2-rt [mH 



= a/2n [mHz] 



Fig. 2. The filled dots represent \k obtained from the 3D 
simulations for the wavenumber k at which E(k) is maxi- 
mum and at the layer where the excitation is maximum in 
the simulation. The results are presented for three simula- 
tions: Procyon (top), the Sun (middle) and aCenB (bot- 
tom). The solid curves represent the Lorentzian form, Eq. 
([U, the dashed curves the Gaussian form Eq. (j4|), and the 
dot dashed curves the exponential form Eq. ([6]). 



from the 3D simulations, and a Lorentzian behaviour is 
never obtained in the present 3D simulations. 

The results of the calculations of CP using both methods 
are presented in Fig.Ufor the six most representative simu- 
lations. In order to remove the large scattering in the direct 
calculations, we perform a running mean over five frequency 
bins. The results of this averaging are shown by dot-dashed 
lines. The choice of five frequency bins is somewhat arbi- 
trary. However we notice that between 2 to 10 frequency 
bins, the maximum and the shape of the spectrum do not 
significantly change. 



Fig. 3. Same as Fig. [2] for k — 10 fc m i n where fc m i n is the 
first non-zero wavenumber of the simulation. 



Comparisons between direct and semi-analytical calcu- 
lations using either x^ or xf a U show systematic differ- 
ences: the excitation rates obtained with the direct calcu- 
lations are systematically lower than those resulting from 
the semi-analytical method. These systematic differences 
are likely due to the too low spatial resolution of the 3D 
simulations which are used here (see Sect. 15.21 below). 

At high frequency, the use of xf instead of Xk results 
in larger 7 for all stars. This arises from the fact that Xk 
spreads slightly more energy at high frequency than Xk 
does (see Fig. [5]). 

The largest difference between the two types of calcu- 
lation (direct versus semi-analytical) is seen in the case of 
Procyon. Indeed, the simulation of Procyon shows a pro- 
nounced depression around v ~ 1.5 mHz. Such a depression 
is not seen in the semi-analytical calculations. The origin 
of this depression has not been clearly identified yet but 
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is perhaps related to some interference between the turbu- 
lence and the acoustic waves which manifests itself in the 
pressure fluctuations in the 3D work integral but is not 
included in the semi-analytical description. 

5.2. Influence of the 3D simulation characteristics 



Lorentzian 
253 2 53x163 
125x125x82 
50x50x82 




Lorentzian 

Gaussian 

253 2 53x163 




[rnHz 



citation rates according to the semi-analytical method as- 
suming either a Lorentzian x.k or a Gaussian x.k- 

As shown in Fig. [5] (top), the excitation rates computed 
according to the direct calculation increase as the spatial 
resolution of the 3D simulation increases. The excitation 
rates computed with the 3D simulations with the two high- 
est spatial resolutions reach approximately the same mean 
amplitude level, indicating that this level of spatial resolu- 
tion is sufficient for the direct calculations. 

We note that as the spatial resolution increases, the 
semi-analytical calculations using a Lorentzian Xk decrease 
by a factor ~ 2 (not shown here). The differences in the 
semi- analytical calculations based on the 253 x 253 x 163 
simulation and the 125 x 125 x 82 simulation are found 
very small, indicating that this level of spatial resolution is 
sufficient for the semi-analytical calculations too. 

Finally, we note that the excitation rates obtained for 
the 50 x 50 x 82 solar simulation (SI) are approximatively 
two times smaller than excitation rates for the 50 x 50 x 82 
solar simulation otherwise used throughout this work (SO 
hereafter). This difference is attributed to the fact that the 
two simulations do not correspond to the same realization. 
Indeed, as a test, we have extended the duration of the 
simulation SI up to 500 min. The full time series has then 
been divided into subsets of equal duration of 100 min and 
p-mode excitation rates have been computed for each sub- 
set. We find that the maximum in the p-modes excitation 
rates fiy) oscillates from a subset to another about a mean 
value. The observed variations are large: the maximum in 
f(v) can be larger (smaller resp.) by ~ 1.5 (0.5 resp.) times 
the maximum in the power spectrum obtained by averag- 
ing the power spectra of all subsets. Hence we find that at 
low spatial resolution, different realizations yield excitation 
rates that are scattered about a mean value at each fre- 
quency. This dispersion is likely to be responsible for the 
factor of two difference between the excitation rate max- 
ima obtained for the two realizations SO and SI. This type 
of dependency of S° — with the starting time of the time 
series and its duration — is expected to be smaller for sim- 
ulations with resolution higher than 50 x 50 x 82, because 
of the larger number of excitation sources there. This will 
be studied in a subsequent work. 



Fig. 5. Top: As in Fig. [4] for solar simulations only. The 
solid line corresponds to the semi-analytical calculations 
based on a Lorentzian Xk and a simulation with a spa- 
tial resolution of 253x253x82. The other lines are run- 
ning means over five frequencies of the direct calculation 
based on solar simulations with different spatial resolution: 
253x253x82 (dashed line), 125x125x82 (dot dashed line) 
and 50x50x82 (dot dot dashed line) . Bottom: The solid 
and dashed lines have the same meaning as in the top fig- 
ure. The dot-dashed line corresponds to the semi-analytical 
calculations based on a Gaussian Xk- 

In order to assess the influence of the spatial resolution 
of the simulation on our results, we have at our disposal 
three other solar 3D simulations, with a grid of 253 x 253 x 
163, 125 x 125 x 82 and 50 x 50 x 82 (hereafter SI), and a 
duration of ~ 42 min, 70 min and 100 min, respectively. 

We have computed the p-modes excitation rates accord- 
ing to the direct method for those three simulations. For 
each of those simulations we have also computed the ex- 



5.3. Eddy-time correlation: Lorentzian versus Gaussian 

As seen in Sect l5.2l above. the characteristics of the simula- 
tions influence the semi-analytical calculations of the mode 
excitation rates (through the input parameters which en- 
ter the semi-analytical calculations and which are taken 
from the 3D simulation). We want to compare the results 
of the semi-analytical calculations using Xk with the semi- 
analytical calculations using Xk ■ It is then necessary to in- 
sert the 3D inputs in these calculations coming from simu- 
lations with the highest quality, here the highest available 
resolution. 

Fig. [5] (bottom) compares semi-analytical calculations 
using a Lorentzian Xk with those using a Gaussian Xk- All 
theses semi-analytical calculations are here based on the en- 
ergy spectrum of the simulation with the spatial resolution 
of 253 x 253 x 163 (see Sect. I5T2]) . 

The average level of the excitation rates calculated 
according to the direct method and with the simulation 
with the highest spatial resolution is in between the semi- 
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analytical calculations based on Lorentzian Xk and those 
based on a Gaussian Xk, nevertheless they are in general 
slightly closer to the semi-analytical calculations based on 
Lorentzian Xk- This result is discussed in Sect. 16.2.11 

5.4. Maximum ofT as a function of L/M 

Fig. [S] shows T max , the maximum in T, as a function of 
L/M for the direct and the semi-analytical calculations. 

The same systematic differences between the direct and 
the semi-analytical calculations as seen in Fig. 2] are of 
course observed here. Note that the differences slightly de- 
crease with increasing values of L/M. 

We have also computed the excitation rate with the 
semi-analytical method using x\ ■ The maximum excita- 
tion rate as evaluated with x\ is systematically larger than 
both the direct calculations and the semi-analytical results 
based on x5? or xjj?. 

In the solar case, JVax is found to be closer to the value 
derived from recent helioseismic data (jBaudin et all [2005) 
when using a Lorentz ian compared to a Gaussian (see also 
Bclk acem et all [20061 B06b hereafter). The 'observed' ex- 
citation rates are derived from the velocity observations V 
as follows: 



y = 2irT v M(h) V 2 



(9) 



where M is the mode mass, V is the mode velocity ampli- 
tude and h is the height above the photosphere where the 
mode mass is evaluated. The mode line width at half max- 
imum in Hz, r„ = rj/ir, (ji is the mode amplitude damping 
rate in s _1 ) is determined observationally in the solar case. 
Using the recent hel ioseismic measurements of V and T v by 
iBaudin et al.l (|2005h and the mode mass co mputed here for 
our s olar model at the height /i=340 km (cf. IBaudin et all 
12505ft . we find ? max , Q = 6.5 ± 0.7 x lO^ergs" 1 . This 
value must be compared with those found with x\ an d 
Xf, namely T^ ax = 4.9 x lO^ergs" 1 
1.2 x 10 22 ergs _1 respectively. 



and CP^ax ( . 



Scaling laws: All sets of calculations can be reasonably 
well fitted with a scaling law of the form y max cx (L/M) s 
where s is a slope which depends on the considered set of 
calculations. Values found for s are summarized in Table [3] 

• For the semi-analytical calculations, we find s = 2.6 
using Xk> s = 3-0 using Xk an d s = 3.1 for the Gaussian 
form. 

The Lorentzian form results in a power law with a 
smaller slope than the Gaussian. This can be understood 
as follows: A Gaussian decreases more rapidly with v than 
a Lorentzian. As the ratio L/M of a main sequence star 
increases, the mode frequencies shift to lower values. Hence 
p-modes of stars with large values of L/M receive relatively 
more acoustic energy when adopting a Gaussian rather 
than a Lorentzian Xk- It is worthwhile to note that even 
though the ratio L/M is the ratio of two global stellar 
quantities, it nevertheless characterizes essentially the stel- 
lar surface layers where the mode excitation is located since 
L/M cx T c %/g. 

• For the set of direct calculations, some scatter exists as 
a consequence of the large statistical fluctuations in CP max 
and a linear regression gives s = 3.4. As expected, this 
value is rather close to that found with the semi-analytical 
calculations using either Xk or Xk- 
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Fig. 6. T max versus L/M where L is the luminosity and 
M is the mass of the ID models associated with the 3D 
simulations. The triangles correspond to the direct calcula- 
tions (labeled as 'DirEx 3D' in the legend), and the other 
symbols correspond to the semi-analytical calculations us- 
ing the three forms of Xk '■ the crosses assume a Gaussian, 
the diamonds an exponential and the squares a Lorentzian, 
respectively. Each set of CP max is fitted by a power law of 
the form (L/M) s where s is the slope of the power law. The 
line-styles correspond to the three semi-analytical cases and 
the direct calculations, as indicated in the lower right cor- 
ner of the plot. 



method 



direct 
semi-analytical 
semi-analytical 
semi-analytical 



3.4 

Gaussian 3.1 1.0 

exponential 3.0 0.9 

Lorentzian 2.6 0.7 



Table 3. Values found for the slopes s (see Sect. I5.4j) and 
sv (see Sect. 15. 5|) . 'Method' is the method considered for 
the calculations of 7. 



5.5. Maximum of the mode amplitudes (V max ) as a function 

of L/M 

The theoretical oscillation velocity amplitudes V can be 
computed according to Eq. §§§ The calculation requires the 
knowledge of the excitation rates, 7, damping rates, 77, and 
mode mass, M. Although it is possible -in principle- to 
compute the convective d ampings from the 3D simulations 
(jNordlund fc Steinll200lh . it is a difficult task which is un- 
der progress. Ho wever , using for instance Gough's Mixing- 
Length Theory (|1976t 119771 G'MLT hereafter), it is pos- 
sible to compute i] and 7 for different stellar models of 
given L, M and deduce F max , the maximum of the mode 
amplitudes, as a function of L/M at the cost of some in- 
consistencies. 

In ISamadi et al.l (|200l[ ). calculations of the damping 
rates r\ based on G'MLT were performed for stellar models 
with different values of L and M. Although these stellar 
models are not the same as those considered here, it is still 
possible, for a crude estimate, to determine the dependency 
of Vmax with L/M. 
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Henc e we proceed as follow s: For each stellar model com- 
puted in I Samadi et alj (|200lD . we derive the values of r\ and 
M at the frequency ^ max at which the maximum amplitude 
is expected. From t he stars for which so l ar-like oscillations 
have been detected. iBedding fc Kieldsenl (|2003l ) have shown 
that this frequency is proportional to the cut-off frequency. 
Hence we determine f max = (y c / ^ Ci ©) J-Wx,© where v c is the 
cut-off frequency of a given model and the symbol refers 
to solar quantities (t'max,© — 3.2 mHz and f Cj © ~ 5.5 mHz). 
We then obtain (?7 max M max ) as a function of L and M. 

On the other hand, in Sect. 15.41 we have established 
J'max as a function of L and M. Then, according to Eq. 
([9]), we can determine V max (L, M) for the different power 
laws of T max . 

We are interested here in the slope (i.e. variation with 
L/M) of Vmax and not its absolute magnitude, therefore we 
scale the theoretical and observed Vmax with a same nor- 
malization value which is taken as the solar v alue Vmax,© = 
33.1 ± 0.9 cms -1 as determined recently by iBaudin et al.l 
(UK). 

We find that Vmax increases as (L/M) sv with different 
values for sv depending on the assumptions for \k- The 
values of sv are summarized in Table [3] and illustrated in 
Fig. El We find sv ~ 0.7 with x\ an d sv ~ 1.0 with Xk ■ 

These scaling laws must be compared with observations 
of a few stars for which solar-like oscillations have been de- 
tected in Doppler velocity. The observed V max are taken 
from Table 1 of HG02, except for 77 Boo, ( Her A, (3 Vir, 
H P 49933 and ji Ara, f or which we use th e V max quoted 
bvlCarrier et alj (l2003h. IMartic et all (l2001l). IMartic et all 
(|2004h . iMosser et al.l (|2005D and iBouchv et alj (120051) re- 
specti vely and e Oph and 77 Ser quoted by iBarban et al.1 
(|200l . 

Fig. ([7]) shows that the observations also indicate a 
monotonic logarithmic increase of V max with L/M despite 
a large dispersion which may at least partly arise from dif- 
ferent origins of the data sets. For the observations we find 
a 'slope' sv ~ 0.7. This is close to the theoretical slope 
obtained when adopting x\ and definitely lower than the 
slopes obtained when adopting Xk or adopted by HG02. 

6. Summary and discussion 

One goal of the present work has been to validate the model 
of stochastic excitation presented in Paper I. The result of 
this test is summarized in Sect. 16.11 A second goal has 
been to study the properties of the turbulent eddy time- 
correlation, Xk, and the importance for the calculation of 
the excitation rates, CP, of the adopted form of Xk- Section 
deals with this subject. 



6.1. Validation of the excitation model 

In order to check the validity of the excitation model, seven 
3D simulations of stars, including the Sun, have been con- 
sidered. For each simulation, we calculated the p-mode ex- 
citation rates, CP, using two methods: the semi-analytical 
excitation model (cf. Sect. 12. 2[ ) that we are testing, and a di- 
rect calculation as detailed in Sect. [2] In the latter method, 
the work performed by the pressure fluctuations on the p- 
modes is calculated directly from the 3D simulations. 

In the semi-analytical method, CP is computed according 
to the excitation model of Paper I. The calculation uses, as 
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Fig. 7. Same as Fig. [5] for Vm a x/Vnax,©, the maximum of 
the mode amplitudes relative to the observed solar value 
(Vmax,© = 33.1 ±0.9cms _1 ). The filled symbols correspond 
to the stars for which solar-like oscillations have been de- 
tected in Doppler velocity. The lines - except the dashed 
line - correspond to the power laws obtained from the pre- 
dicted scaling laws for CP max (Fig. [5]) and estimated values 
of the damping rates f/ max (see text for details) . Results for 
two different eddy time-correlation functions, Xk, are pre- 
sented: Lorentzian (solid line) and Gaussian (dot-dashed 
line) functions. For comparison the dashed line shows the 
result by HG02. Values of the slope sv are given on the plot 
and in Table 3. 



input, information from the 3D simulations as for instance 
the eddy time-correlation (xk) and the kinetic energy spec- 
tra (E(k)). However although Xk has been computed for 
each simulation, in practice for simplifying the problem 
of implementation as well as for comparison purpose with 
Paper III, we chose to represent the v variation of x.k with 
simple analytical functions. It is found that the z/-variation 
of Xk in the present simulations lies loosely between that of 
an exponential and a Gaussian. We then perform the val- 
idation test of the excitation model using those two forms 
of Xk- 

We find that using either x*k 01 xf m the semi-analytical 
calculations of CP results in systematically higher excita- 
tion rates than those obtained with direct 3D calculations. 
These systematic differences are attributed to the low spa- 
tial resolution of our present set of simulations. Indeed we 
have shown here that using solar simulations with different 
spatial resolutions, the resulting excitation rates increase 
with increasing spatial resolution. 

We have next investigated the dependence of CP max with 
L/M (See Fig. 0, where L and M are the stellar luminos- 
ity and mass respectively. As in pr evious works based o n 
a purely theoretical approach (e.g. ISamadi et al.1 . l2003al ). 
we find that CP max scales approximatively as (L/M) s where 
s is the slope of the scaling law: we find s = 3.4 with the 
direct calculations and s = 3.2 and s = 3.1 with the semi- 
analytical calculations using Xk an d Xk respectively. This 
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indicates a general agreement between the scaling proper- 
ties of both types of calculations, which validates to some 
extent the adopted excitation model across the domain of 
the HR diagram studied here. 

For the sake of simplicity, only simple analytical forms 
for Xk have been investigated here. We expect that the 
use of more sophisticated forms for Xk would reduce the 
dispersion between the analytical and direct calculations, 
but would not affect the conclusions of the present paper. 

6.2. The eddy time-correlation spectra, Xk 

The slope s of the scaling law for y max , is found to de- 
pend significantly on the adopted analytical form for Xk- 
The semi-analytical calculations using the Lorentzian form 
for Xk results in a significantly smaller slope s than those 
based on the Gaussian or the exponential or from direct 
calculations (see Table [3]). 

Except for the Sun, independent and accurate enough 
constraints on both the mode damping rates and the mode 
excitation rates are not yet available. We are then left to 
perform comparison between predicted and observed mode 
amplitudes. Unfortunately, obtaining tight constraints on 
Xk using comparison between predicted and observed mode 
amplitudes is hampered by large uncertainties in the theo- 
retical estimates of the damping rates. It is therefore cur- 
rently difficult to derive the excitation rates D 3 for the 
few stars f or which solar - like o scillations have been de- 
tected (see ISamadi et al.l. l2004f). The fu t ure s pace mis- 
sion CORQT (|Baglin &; The Corot Team! Il998t ) will pro- 
vide high-quality data on seismic observations. Indeed the 
COROT mission will be the first mission that will provide 
both high precision mode amplitudes and line-widths for 
stars other than the Sun. It will then be possible to use 
the observed damping rates and to derive the excitation 
rate 7 free of the uncertainties associated with a theoreti- 
cal computation of damping rates. In particular, it will be 
possible to determine y max as a function of L and M from 
the observed stars. Such observations will provide valuable 
constraints for our models for Xk- 

We can, nevertheless, already give some arguments be- 
low in favor of the Lorentzian being the correct description 
for Xk- 

6.2.1. Solar case 

In the 3D simulations studied here, including that of 
the Sun, the inferred v dependency of Xk is far from a 
Lorentzian, in contrast to that found with the solar 3D 
simulation investigated in Paper III. However, by investi- 
gating solar simulations with different resolutions, we find 
that, as the spatial resolution increases, Xfc tends towards 
a Lorentzian ^-dependency. This explanation is likely to 
stand for non-solar simulations too, but has not yet been 
confirmed (work in progress). 

Furthermore, as shown in Fig. O bottom, the direct 
calculations obtained with the simulation with the high- 
est spatial resolution available is slightly closer to the semi- 
analytical calculations using the Lorentzian form than those 
using the Gaussian one. 

Independently of the resolution (if large enough of 
course), a Lorentzian Xk predicts larger values for y max 
than a Gaussian or an exponential do. In particular in the 



solar case, the semi-analytical calculation using x\ results 
in a JVax clo s er to the helioseismic constraints derived by 
iBaudin et al.l (|2005j ) compared to using Xk or Xk- This lat- 
ter result is in agreement with that of Paper III. 

Part of the remaining discrepancies with the helio- 
seismic constraints are attributed to the a dopted clo- 
sure model according to lBelkacem et all (|2006l B06b here- 
after). Indeed, theoretical models of stochastic excitation 
adopt the quasi- normal approximation (QNA). As shown 
in B06b, the skew introduced by the QNA result in a 
under-estimation of the solar p mode excitation rates. 
W hen the so-called closu re model with plumes proposed 
bv lBelkacem et all (|2006al ) is adopted, new semi-theoretical 
calculations fit ra ther well the recen t helioseismic con- 
straints derived bv lBaudin et al.l (|2005l . see B06b). 

6.2.2. Vmax as a function of L/M 

Consequences of the predicted power laws for CP max have 
also been crudely investigated here for the expected value 
of Vmax, the maximum value of the mode velocity (Fig. [7j). 
Calculations of Vmax from y max require the knowledge of 
the mode damping rates, r/, which cannot be fully deter- 
mined from the simulations. We are then led to use theo- 
retical calculations of the damping rates. We consider here 
those perfor med by ISamadi et al.l (|200lh which are based 
on Gough's (|l976l : Il977h non-local and time-dependent for- 
mulation of convection. From those values of rj and the 
different power laws for y max expected values of Vmax are 
obtained. 

We find, as in iHoudek fc Goughl (|2002ft (HG02), that 
Vmax scales as (L/M) sv . Calculations by HG02 result in 
sv ~ 1.5. Our semi-analytical calculations of Tmax based on 
a Gaussian Xk result in a slightly smaller slope (sv ~ 1.0). 
On the other hand, using a Lorentzian Xk results in a slope 
sv ~ 0.7 which is closer to that derived from the few stars 
for which oscillation amplitudes have been measured. 

From this result, we conclude that the problem of the 
over-estimation of the amplitudes of the solar-like oscil- 
lating stars more luminous than the Sun is related to the 
choice of the model for Xk- Indee d, previous theoretical cal- 
culations bv lHoudek et al.l (|l999h are based on the assump- 
tion of a Gaussian Xk ■ As shown here, the Gaussian assump- 
tion results in a la rger slope sv than the Lorentzian Xk ■ This 
is the reason whv lHoudek et all (|l999f ) over-estimate Vmax 
for L/M > Lq/Mq. 

On the other hand, if one assumes Xk = x\> a scaling fac- 
tor is no longer required to reproduce y max for the solar 
p-modes. Moreover, as a consequence of the smaller slope, 
sv, resulting from a Lorentzian Xkj the predicted ampli- 
tudes for other stars match the observations better. 

This result further indicates that a Lorentzian is the 
better choice for Xk, as was also concluded in Paper III. 

Departures of the theoretical curve from the observed 
points in Fig. [7] can be attributed to several causes which 
remain to be investigated: 

1) A major uncertainty comes from the computed damping 
rates as no accurate enough observations are available 
yet to validate them. As V results from the balance be- 
tween y and 77, the slope sv can also depend on the 
variation of r\ with L/M . Thus, the large differences in 
sv between the seismic observations and the calculations 
based on Xk can also he, a priori, attributed to an in- 
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correct evaluation of the damping rates. However ?7 max 
- the value of the damping rate at the frequency ^ max at 
which the maximum amplitude is expected - does not 
follow a clear scaling law with L/M. We have looked 
at the ?7 max variation in our set of G'MLT models and 
found no clear dependence of 77 max on L/M but rather 
a dispersion. 

2) The observed stars in Fig. [7] have somewhat differ- 
ent chemical compositions; this can cause some scat- 
ter in the relation V meix -L/M which has not been taken 
into account here. All the simulations investigated in 
the present work employ a solar metal abundance. The 
metallicity has a direct impact on the opacity and the 
EOS. Both in turn affect the internal structure and are 
also decisive for the transition from convection to ra- 
diation in the photosphere and therefore determine the 
structure of the super-adiabatic region. Hence, the prop- 
erties of the super-adiabatic region, relevant for the ex- 
citation rates, differ for stars located at the same po- 
sition in the HR diagram (e.g., same T e g and same g) 
but with different metal abundances. Consequently the 
excitation of p-modes for such stars probably differ, al- 
though it remains to be seen to what extent. A differen- 
tial investigation of the metallicity effect is planned for 
the future. 

6.3. Relative contribution of the turbulent pressure 

Another issue concerns the relative contribution of the 
turbulent pressure. The excitation of solar-like oscilla- 
tions is generally attributed to the turbulent pressure (i.e. 
Reynolds stress) and the entropy fluctuations (i.e. non- 
adiabatic gas pressure fluctuations) and occurs in the super- 
adiabatic region where those two terms are the largest. In 
Paper III, it was found that the two driving sources are of 
the sam e order of magnitude, in c ontradiction with the re- 
sults by IStein fc Nordlundl (|200ll ) who found - based on 
their 3D numerical simulations of convection - that the 
turbulent pressure is the dominant contribution to the ex- 
citation. The discrepancy is removed here as we used a cor- 
rected version of the formulation of the contribution of the 
Reynolds stress of Paper I (see Eq. ©), leading to a larger 
contribution from the Reynolds stress. 

For the Sun, assuming \\ (Xk resp.), we now find that 
the Reynolds stress contribution is 5 times (3 times resp.) 
larger than that due to the entropy fluctuations (non- 
adiabatic gas-pressure fluctuations). Hence, the Reynolds 
stress is indeed the domin ant source of excita ti on in agree- 
ment with the results of IStein fc Nordlundl (|2001h . The 
best agreement with the latter results is obtained with a 
Lorentzian \k- 

However, we find that the relative contribution from 
Reynolds stresses decreases rapidly with (L/M). For in- 
stance, in the simulation of Procyon, the Reynolds stress 
represents only ~ 30 % of the total excitation rate. 

From that, we conclude that the excitation by entropy 
fluctuations cannot be neglected, especially for stars more 
luminous than the Sun. 
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Appendix A: Calculation of the non adiabatic 
pressure fluctuations 

The adiabatic variation of the gas pressure does not con- 
tribute to the A(PdV) work over an oscillation period as it 
is in phase with the volume (or density) variation. In prac- 
tice, however, it is beneficial for the accuracy of the com- 
putation of excitation to subtract the adiabatic part of the 
gas pressure fluctuation, since it reduces the coherent part. 
That part gives zero contribution only in the limit of infinite 
time, or for an exact integer number of periods. However, in 
practice, it gives rise to a random (or noisy) contribution. 
Indeed, as we deal with a lot of different modes it is hard to 
find a time-interval which is an integer number of periods 
of each and all of the modes at the same time. 

The Lagrangian variations of gas pressure, AP gas must 
satisfy 

AP gas = Ap + ^AS (A.l) 

p ob 

where -P gas , p and S are the gas pressure the density and the 
entropy respectively and where the operator A represents 
the pseudo Lagrangian fluctuations of a given quantity. The 
concept of pseudo Lagrangian fluctuations is introduced in 
Nordlund & Stein (2001). Accordingly we derive the non- 
adiabatic gas pressure fluctuations as: 

AP gas , nad (r, t) = AP gas - c 2 s Ap (A.2) 

where c 2 s = TiP glis /p is the sound speed. 

However, what we want to subtract off from AP gas is 
that part of the pressure variation that is due to adiabatic 
compression and expansion due to the particular radial 
wave modes (i.e. the low amplitude perturbation of p(r) 
on top of the possibly large variations horizontally of p(r) 
that p(r) is an average of). 

To find the nonadiabatic pressure fluctuations, we start 
with calculations of horizontal averages of the primary 
quantities, P gas , Pturb, P and c 2 . We convert these aver- 
ages to the pseudo-Lagrangian frame of reference, in which 
the net mass flux vanishes. We then compute fluctuations of 
the resulting quantities with respect to time, i.e., subtract 
their time averages: 

AP gas (Pgas)h (-fgas)h,t 
APturb = (Purb)h — (Pturb) h,t 
A/? = (p)h - (P)h,t 

Here, ()h refers to horizontal average and ()h,t refers to con- 
sequent time average performed on a horizontally averaged 
quantity. Finally, the non-adiabatic fluctuations of the total 
pressure (that is gas + turbulent pressure) are: 

AP nac [ — AP gasna d ~\~ APt ur b 

= AP-( C 2) h , t Ap (A.3) 
where AP = AP gas + AP turb . 
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Fig. 4. Excitation rates, CP, are presented as functions of mode frequency for six of the seven convection simulations 
listed in Tables [1] and [51 Each triangle corresponds to a single evaluation of the 3D work integral estimated for a given 
eigenfrequency according to Eq. ([T]) . The dot-dashed lines correspond to a running mean of the triangle symbols performed 
over five frequencies. The solid and dashed lines correspond to the excitation rates calculated with the semi-analytical 
method and using the Gaussian and the exponential forms of \k, respectively. All results shown are obtained as the 
sum of contributions from the two sources of excitation: excitation by the turbulent pressure and excitation by the 
non-adiabatic gas pressure. 



